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1.  Introduction 


This  document  serves  as  the  final  written  report  for  agreement  DAAD17-01-P-1238  between  the 
U.S.  Army  Research  Laboratory  (ARL)  and  the  Alabama  A&M  University  Research  Institute.  It 
documents  adherence  to  the  following. 

1.1  Statement  of  Work 

Development  of  the  Message  Passing  Interface  -based  Visual  extensible  Molecular  Dynamics 
(VXMD)  compatible  modules  for  the  creation,  maintenance,  and  migration  of  parallel  atomic 
bond  lists  and  the  force  modules  that  use  the  bond  lists. 

The  principal  investigator  is  Dr.  Andrew  Scott. 

The  work  under  this  agreement  will  consist  of  extending  and  improving,  and  parallelizing  the 
serial  bond  list  software  written  by  Dr.  Betsy  Rice.  The  extended  software  will  have  two 
components,  one  for  handling  the  atomic  bond  lists  within  the  domain  decomposition  used  by 
VXMD,  and  the  second  will  be  the  force  and  integration  routines  that  use  these  lists.  These  lists 
require  the  development  of  a  mapping  scheme  that  will  transform  “skewed”  cells  into  rectangular 
for  two-dimensional  (2-D)  and  three-dimensional  (3-D)  cases. 

The  specific  deliverables  for  this  task  will  be  a  robust  set  of  modules  that  allow  the 
implementation  of  a  MD  simulator  assuming  the  NPT  ensemble. 

All  results  and  codes  will  be  shared  with  ARL.  The  final  document  will  describe  algorithm  and 
code  development,  as  well  as  specifics  for  interfacing  with  existing  software  under  CCM-5.  Any 
follow-on  computational  effort  will  focus  on  extending  the  VXMD.  Other  issues  will  be 
addressed  as  they  arise. 

1.2  Period  of  Performance 

This  report  describes  work  that  was  conducted  from  1  October  2001-3 1  January  2002. 
Preliminary  results  were  provided  to  ARL  in  the  November  2001  timeframe. 

Specific  deliverables: 

•  Routines  that  transform  non-rectangular  prism  simulation  cells  to  rectangular  prism  cells, 

•  Within  shape-changing  NPT  MD  simulation,  coordinate  transforms  such  that  optimal 
rectangular  prism  simulation  shape  is  maintained,  incorporated  into  the  MPI  version  of  the 
MD  code,  and 
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•  For  post-simulation  analysis,  routines  for  inverse  transformation  of  simulation  cell  from 
rectangular  prism  shape. 

This  report  provides  an  overview  of  the  algorithms  developed  and  the  software  produced,  and 
how  they  were  incorporated  in  the  existing  MD  software.  The  actual  source  code  was  submitted 
to  Dr.  Betsy  Rice  at  ARL  on  30  January  2002.  The  new  code  that  was  generated,  and  the 
existing  Novice-MD  routines  that  were  modified,  are  detailed  in  section  3.* 


2.  Algorithm  Overview 


The  following  describes  a  methodology  for  converting  the  global  Cartesian  lattice  system  to  a 

“wrapped”  local  cell  coordinate  system. 

Given: 

•  Lattice  vectors  [a,b,c]  in  global  Cartesian  fonnat. 

Determine: 

•  Coordinate  transfonnation  matrix,  A,  to  map  the  atoms  and  lattice  vectors  into  the  local 
simulation  environment  for  the  initial  lattice  and  at  each  time  step. 

•  “Wrap”  the  resultant  local  lattice  vectors  into  rectangular  (right  angled)  unit  cells  for 
solution  by  the  MD  code  for  the  initial  lattice  and  at  each  time  step. 

•  Track  the  global  transformation,  Agi0bai,  to  provide  conversion  back  to  the  initial  global 
coordinates  at  the  end  of  the  simulation,  and/or  at  intervals  throughout. 

Products: 

•  The  first  lattice  vector,  a,  must  be  oriented  in  the  local  x-direction. 

•  The  next  lattice  vector,  b  (considering  only  the  components  orthogonal  to  the  first  lattice 
vector)  must  be  mapped  into  the  local  y-direction. 

•  The  third  lattice  vector,  c  (considering  only  the  components  orthogonal  to  the  previous 
lattice  vectors)  must  be  mapped  into  the  local  z-direction. 

Procedure: 

1 .  Read  in  lattice  vectors  [a,b,c], 

2.  Find  a  set  of  global  orthogonal  vectors  [la,lb,lc]  from  [ a,b,c ]. 


*For  further  information  regarding  dissemination  of  the  software,  contact  Dr.  Betsy  Rice,  USARL,  ATTN:  AMSRD-WM- 
BD,  Aberdeen  Proving  Ground,  MD  21005-5066,  <betsyr@arl.army.mil>. 
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(a)  Denote  first  cell  vector  11. 

(1)  Global  coordinates  of  11  match  the  global  coordinates  of  the  first  lattice  vector,  a. 

(2)  Determine  the  global  unit  vector  of  11,  llunit. 

(b)  Determine  the  next  cell  vector  (orthogonal  components  to  11)  from  the  lattice  vector,  b. 

(1)  Find  the  magnitude  of  the  components  parallel  to  //;  denote  as  b_shadow,  e.g.,  b_shadow 
=  dot _product(b,ll  unit) . 

(2)  Find  the  components  of  b  that  are  parallel  to  11;  denote  as  b _par,  e.g.,  b _par 
=  b_shadow*I  limit. 

(3)  Find  the  components  of  b  that  are  perpendicular  to  11;  denote  as  b_orth,  e.g.,  b  orth  =  b  - 
b  _par=12. 

(4)  Note  that  the  global  coordinates  of  12  represent  only  the  components  of  the  original  lattice 
vector,  b,  that  are  orthogonal  to  11 . 

(5)  Determine  the  global  unit  vector  of  12, 12un.it. 

(c)  Determine  13,  from  the  remaining  lattice  vector,  c. 

(1)  Find  the  components  of  c  that  are  parallel  to  11,  as  above.  Denote  c _parl. 

(2)  Similarly,  find  the  components  of  c  that  are  parallel  to  12  as  above.  Denote  c _par2. 

(3)  Determine  the  components  of  c  that  are  perpendicular  to  both  11  and  12.  Denote  13  =  c  - 
c _parl  -  c _par2. 

(d)  [11,12,13]  are  orthogonal  cell  vectors  in  global  coordinates  formed  by  projection  of  the 
initial  lattice  vectors. 

3.  Find  a  set  of  local  cell  vectors  [11  iocai,  l21ocal,  131ocal]  from  the  global  vectors  [11,  12,  13]. 

(a)  Find  the  magnitudes  (lengths)  of  the  vectors  \ll\,  \12\,  |/3|. 

(b)  Map  the  magnitude  of  the  11  vector  to  the  local  x-direction. 

(c)  Map  the  magnitude  of  the  12  vector  to  the  local  y-direction. 

(d)  Map  the  magnitude  of  the  13  vector  to  the  local  2-direction. 

(e)  For  example,  the  mapping  from  the  global  [11,12,13]  to  the  local  system 

[U local, 12 local,  13 local]  takes  the  following  form: 
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4.  Find  the  coordinate  transformation  matrix,  A. 

(a)  Denote  the  set  of  global  vectors  [11,12,13]  as  the  matrix  L  and  the  set  of  local  vectors 

[l  1  local]2 lOCal,l3 local]  &S  the  matrix  Llocal • 

(b)  The  mapping  function  takes  the  form  of  the  equation,  AL  =  Liocai,  where  A  is  an  unknown 
coordinate  transformation  matrix  providing  the  mapping  from  the  global  to  the  local 
system. 

5 .  Solve  for  A,  with  the  following: 

A  =  Lheairl.  (1) 

6.  Convert  global  coordinates  to  local. 

(a)  Determine  the  local  lattice  vectors,  with  A*[a,b,c]  =  [aiocai,biocai,cioca] . 

(b)  Determine  the  local  atom  positions  with  A  *rgi0bai  =  r local¬ 
ise)  Etc.... 

7.  Determine  the  form  of  the  “wrapped”  rectangular  unit  cell. 

(a)  The  global  matrix  L  and  the  local  matrix  L/ocai  meet  the  criteria  for  the  “wrapped” 
rectangular  unit  cell  vectors  since  they  are  orthogonal  and  created  from  projections  of  the 
initial  lattice  vectors. 

(b)  This  can  be  verified  by  checking  that  the  volume  of  the  “wrapped”  cell  is  equal  to  that  of 
the  original  lattice  cell. 

(c)  Atom  coordinates  can  be  wrapped  to  fit  into  the  rectangular  unit  cell  by  utilizing  the 
assumption  of  periodic  conditions  in  the  overall  lattice  structure. 

8.  Update  the  historical  transformation  matrix,  which  can  be  used  to  convert  back  to  the 
original  or  “global”  coordinate  system  when  required. 

(a)  If  t=t(0),  then  Ahlstory  =  A. 

(b)  Else  A  history  A  *  A  history- 

In  the  case  of  NPT/ShapeChange,  there  is  a  strain  tensor  acting  on  the  cell  vectors  at  each 
iterative  step  in  the  simulation.  The  actual  lattice  vectors  need  to  be  updated  with  this 
incremental  strain  in  order  to  identify  specific  boundary  conditions.  The  lattice  vectors  in  local 
coordinates  are  updated  as  follows: 
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•  Determine  the  incremental  strain  in  the  simulation  cell,  dStrainTensor,  by  multiplying  the 
strained  local  cell  vectors,  L  local  strained,  by  the  inverse  of  the  original  local  cell 
vectors,  L_localn, 

dStrainTensor  =  (  L_local_strained)*(L_Localn  .  (2) 

•  Detennine  the  new  (strained)  actual  lattice  vectors  in  local  coordinates  by  multiplying  the 
current  lattice  vectors  by  the  newly  acquired  dStrainTensor. 

latticeVectorsn+i  =  dStrainTensor*  latticeVectorsn.  (3) 

In  the  case  of  RungeKuttaNPT,  a  shape  change  is  required  at  each  of  the  four  steps  in  the 
integration  method.  Therefore,  the  transfonnation  is  applied  to  each  of  the  relevant  variables  in 
each  step  and  any  corresponding  variables  obtained  in  previous  steps.  This  is  required  since  the 
Runge-Kutta  method  adds  terms  from  each  of  the  four  steps,  so  they  must  be  represented  in  the 
same  coordinate  system  to  provide  meaningful  results. 

In  the  case  of  LeapFrogVerletNPT,  a  shape  change  is  required  only  once  in  the  integration 
method.  It  applies  to  any  active  variables,  as  well  as  historical  values  that  are  saved  for  use  in 
subsequent  iterations. 


3.  Code  Overview 


The  following  describes  new  code  that  was  written  and  outlines  the  changes  that  were  made  in 
the  existing  software  to  accommodate  the  new  functionality.  The  description  is  listed  by 
individual  source  code  file  and  routine  name  as  appropriate.  The  Fortran  source  files  are  denoted 
with  the  .f90  suffix.  Individual  routines  and/or  data  objects  that  are  being  described  are  listed  in 
italics.  Items  in  single  quotes  refer  to  variables  and  routines  within  the  Novice-MD  software. 
Items  in  double  quotes  refer  to  system  variables  and  routines.  The  double  colon  convention  is 
characteristic  of  software  specification.  As  previously  mentioned,  contact  Dr.  Betsy  Rice  for 
further  infonnation  regarding  dissemination  of  the  software. 

accumulateMod.f90 

•  IncrementalOutput: :  rotate  incremental  stress  tensor  to  global  coordinates.  Call 
TncrementAccumulator’.  Rotate  ‘stress  tensor’  back  to  local  coordinates. 

boundaryConditionsMod.f90 

•  BoundaiyConditionsUpdate::  if  ‘ShapeChange’,  use  the  ‘current_lattice_  vectors’  from  the 
transformCell.  Otherwise  use  the  simulationCell  lattice  vectors. 

brennerMod.f90 
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•  Brenner  Initialize: :  Fixed  bug,  initialize  ‘dimension’  variable. 
controlDescMod.f90 

•  ControlDescCreate: :  changed  “rand”  and  “srand”  functions  to  “ranf  ’  and  “seed”  functions 
for  use  on  the  CRAY  platform.  SGI  requires  the  original  file. 

latticeMod.fDO 

•  SimulationCellRestart: :  if  ‘ShapeChange’  then  read  the  transformCell  from  the  restart  file 
with  routine  ‘SimulationTransfonnRead’. 

•  LatticeNewRead: :  if  the  root  process,  rotate  the  atoms  into  the  local  system  using 
‘UtilityRotateCoordinates’,  then  wrap  them  into  the  local  simulationCell  using 
‘Utility  WrapCoordinates’. 

•  LatticeWrite: :  if  ShapeChange’,  output  the  transformCell  to  the  restart  file  using 
‘SimulationTransfrom  Write’. 

•  LatticeAtomsPrint: :  if  the  root  process,  then  rotate  the  atom  coordinates  into  global  space, 
write  the  data  out,  then  rotate  the  atom  coordinates  back  into  local  space  using 
‘UtilityRotateCoordinates’. 

•  Fee  and  cubic::  changed  ‘rand’  to  ‘ranf  for  use  on  the  CRAY;  reverse  for  SGI. 
moverMod.f90 

•  LeapFrogVerletNPT: :  is  ‘ShapeChange’,  incorporate  transformCell  to  allow  shape  change 
of  the  current  cell  and  subsequent  cell  rectangularization  for  each  iteration,  using 
‘SimulationCellUpdate’.  Rotate  coordinate  variables;  x,  vx,  dvx,  void  into  new  local 
coordinates  using  ‘UtilityRotateCoordinates’.  Rotate  the  strain  matrix,  eta,  into  the  new 
system. 

noviceMD.f90  (main  program) 

•  Create  the  ‘transformCell’  from  the  initial  ‘simulationCell’  input  using 
‘  S  imulationCellT  ransfonn  ’ . 

•  Set  ShapeChangeFlag  =  1  if  ‘ShapeChange’. 

•  Added  ‘DEBUG’  statements  to  print  transformCell. 
rungeKuttaMod.f90 

•  RungeKuttaNPT: :  if ‘ShapeChange’,  incorporate  ‘transformCell’  to  allow  shape  change  of 
unit  cell  and  subsequent  cell  rectangularization  for  each  of  the  four  steps  in  the  R-K 
integration  using  ‘SimulationCellUpdate’.  Rotate  coordinate  variables;  x,  vx,  dvx,  kO,  kl, 
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10,  11,  q,  and  qv  for  each  step  using  ‘UtilityRotateCoordinates’.  Rotate  the  matrices;  eta, 
eO,  el,  hO,  hi  into  new  local  coordinate  system. 

simulationCellMod.f90 

•  SimulationCellT : :  added  a  pointer,  ‘transfonnCell’  of  type  ‘transfonnCellT’  to  the  data 
structure  for  ready  access  to  the  transformation  data. 

•  SimulationCeUCreate::  initialize  ‘transfonnCell’  pointer  to  NULL. 

•  SimulationCellUPdate: :  added  this  routine  to  update  the  transform  and  simulation  cell 
structures  during  NPT  ‘ShapeChange’  procedure. 

•  SimulationCell  Transform: :  added  this  routine  to  transform  the  initial  simulationCell 
lattice  vectors  to  the  rectangularized  cell  for  simulation. 

•  SimulationCellFormattedWrite: :  if  ‘ShapeChange’  then  transform  current  lattice  vectors  to 
global  coordinates.  Output  global  lattice  vectors.  Transform  lattice  vectors  back  into  local 
coordinate  system. 

•  SimulationTransformWrite: :  added  this  routine  to  write  the  transfonnCell  structure  to  the 
‘RESTART’  file,  in  the  case  of ‘ShapeChange’. 

•  SimulationTransformRead: :  added  this  routine  to  read  the  transfonnCell  structure  from  the 
‘RESTART’  file. 

•  SimulationCellDestroy::  dispose  of ‘transfonnCell’  structure. 

•  CalculateNType: :  changed  the  function  for  random  number  generation  to  RANF  for 
CRAY  port,  and  back  to  RAND  for  SGI. 

transformCellMod.f90 

•  TransformCellT: :  added  the  data  structure  containing  the  components  needed  to  transform 
the  ‘simulationCell’  to  rectangularized  cells.  Includes  rotational  history  data  and 
‘strainTensor’  history  infonnation. 

•  TransformCellCreate: :  new  routine  to  allocate  memory  and  initialize  the  ‘transfonnCell’ 
data  structure. 

•  TransformCellUpdate: :  new  routine  to  update  the  ‘transfonnCell’  structure  given  a  new 
lattice  vector  configuration. 

•  TransformCellProcess: :  Performs  the  actual  calculations  to  find  the  rotated,  wrapped  cell 
and  the  current  transform  matrix  A. 

•  TransformCellPrint: :  new  routine  to  release  the  memory  associated  with  the 
TransfonnCell  structure. 
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utilityMod.f90 

•  Utility RealPrint: :  new  routine  to  output  a  formatted  array  of  real  numbers.  Typically 
some  sort  of  coordinate  information. 

•  Utility  WrapCoordinates: :  new  routine  to  “wrap”  an  array  of  atom  coordinates  to  fit  inside 
a  given  simulation  cell. 

•  UtilityRotateCoordinates::  new  routine  to  rotate  an  array  of  atom  coordinates  with  a  given 
rotation  transform  matrix,  A. 

•  Utility  Inverse: :  new  routine  to  “wrap”  an  array  of  atom  coordinates  to  fit  inside  a  given 
simulation  cell. 

•  UtilityMag: :  new  function  to  calculate  the  magnitude  of  a  given  vector. 

•  UtilityUnit: :  new  routine  to  determine  the  unit  vector  associated  with  a  given  vector 

•  UtilityPar: :  new  routine  to  detennine  the  parallel  components  of  one  given  vector  to 
another. 

•  UtilityPerp: :  new  routine  to  detennine  the  perpendicular  components  of  one  given  vector 
to  another. 

•  UtilityCross_3d::  new  routine  to  determine  the  vector  cross  product  of  two  3-D  vectors. 

•  UtilityGauss::  Public  domain  f77  LINPACK  LU  decomposition  routine,  which  was 
modified  with  f90  syntax  and  customized  for  use  with  this  project. 

•  UtilityGedi : :  Public  domain  f77  LINPACK  routine,  which  inverts  the  LU  decomposed 
matrix  resultant  of  UtilityGauss  and  calculates  its  detenninant  (used  for  cell  volume 
calculation).  It  was  modified  with  f90  syntax  for  this  project. 

•  UtilityScal: :  Public  domain  f77  BLAS  routine  that  scales  a  given  vector  by  a  constant.  It 
was  modified  with  f90  syntax  for  this  project. 

•  UtilitySaxpy: :  Public  domain  f77  BLAS  routine  that  calculates  a  constant  times  a  vector 
plus  a  vector.  It  was  modified  with  190  syntax  for  this  work. 

•  UtilitySwap::  Public  domain  f77  BLAS  routine  that  interchanges  two  vectors.  It  was 
modified  with  f90  syntax  for  this  work. 

•  Utilitylsamax: :  Public  domain  f77  BLAS  routine  that  returns  the  index  number  of  the 
element  having  the  maximum  value  in  a  given  real  array.  It  was  modified  with  f90  for  use 
in  this  project. 

initialV  elocitiesMod.f90 
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•  Gauss::  changed  random  number  generator  to  “RANF”  for  use  on  the  CRAY  and  back  to 
“RAND”  for  the  SGI. 

mallopt.c 

•  MALLOPTDEBUG/MALLOPTNODEBUG::  added  an  extern  “C”  declaration  for  use  on  the 
CRAY.  Changed  the  ‘MDEBUG’  keyword  to  ‘M_TRA.CE’  on  the  CRAY  and  back  for 
SGI.  The  keywords  differ  between  the  platforms. 

putAttributeMod.f90,  conserveMod.f90,  ewaldMod.f90,  potentialMod.f90,  radial. f90, 

radialMod.f90,  rungeKuttaMod.f90,  workMod.190 

•  Changed  the  “double  precision”  variable  type  to  “real  (double)”  for  consistency  and  porting 
to  the  CRAY. 

makefile  (/Cells  directory) 

•  Modified  the  ‘makefile’  to  include  the  new  source  code  files,  transformCellMod.f90  and 
utilityMod.f90  for  inclusion  in  the  compilation  process. 

•  Updated  the  dependency  map  to  ensure  appropriate  ordering. 


4.  Verification 


The  methodology  and  concomitant  code  were  tested  and  verified  using  one-dimensional,  2-D, 
and  3-D  cell  vectors.  The  matrix  inversion  and  vector  routines  are  dimensionally  independent, 
with  the  exception  of  “UtilityCross-3d”,  which  requires  3-D  vector  input. 

The  updated  program  was  compiled  and  tested  on  CRAY  SV1  and  SGI  0rigin2000  platforms. 
The  preliminary  tests  indicate  that  the  desired  ‘NPT/ShapeChange’  functionality  is  observed. 


5.  Deliverables 


The  source  code  for  both  CRAY  and  SGI  versions  were  electronically  submitted  to  Dr.  Rice  on 
30  January  2002.  This  report  satisfies  the  written  documentation  necessary  for  completion  of  the 
project.  As  mentioned  previously,  Dr.  Rice  may  be  contacted  for  information  regarding 
dissemination  of  the  software. 
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6.  Future  Work 


Subject  to  the  availability  of  continued  funds,  future  work  will  include  updating  the  cell’s 
neighbor  lists  and  incorporating  parallel  force  communications  to  the  MD  simulations,  or  any 
other  work  deemed  necessary  by  Dr.  Rice. 


7.  Conclusions 


This  report  has  outlined  the  conformance  and  production  of  the  project  deliverables,  specifically: 

•  Routines  that  transfonn  non-rectangular  prism  simulation  cells  to  rectangular  prism  cells 
were  created  and  tested. 

•  Within  shape-changing  NPT  MD  simulation,  coordinate  transformations  to  maintain  the 
optimal  rectangular  prism  simulation  shape  were  developed  and  introduced  into  the  MPI 
version  of  MD  code. 

•  For  post-simulation  analysis,  routines  for  inverse  transformation  of  simulation  cell  from 
rectangular  prism  shape  were  created. 

In  addition,  the  MD  code  was  ported  and  tested  on  the  Cray  SV-1  vector  platform. 
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NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


ABERDEEN  PROVING  GROUND 
1  DIR  USARL 

AMSRD  ARL  Cl  OK  TP  (BLDG  4600) 


1  US  ARMY  RSRCH  DEV  & 
ENGRG  CMD 
SYSTEMS  OF  SYSTEMS 
INTEGRATION 
AMSRD  SS  T 
6000  6TH  ST  STE  100 
FORT  BELVOIR  VA  22060-5608 

1  INST  FOR  ADVNCD  TCHNLGY 
THE  UNIV  OF  TEXAS 
AT  AUSTIN 

3925  W  BRAKER  LN  STE  400 
AUSTIN  TX  78759-5316 

1  US  MILITARY  ACADEMY 

MATH  SCI  CTR  EXCELLENCE 
MADN  MATH 
THAYER  HALL 
WEST  POINT  NY  10996-1786 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
IMNE  AD  IM  DR 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 

3  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRD  ARL  Cl  OK  TL 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 

3  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRD  ARL  CS  IS  T 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 


1  DEFENSE  TECHNICAL 
(PDF  INFORMATION  CTR 
ONLY)  DTICOCA 

8725  JOHN  J  KINGMAN  RD 
STE  0944 

FORT  BELVOIR  VA  22060-6218 
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NO.  OF 

COPIES  ORGANIZATION 

1  ALABAMA  A&M 

DEPT  ELEC  ENGNRNG 
A  SCOTT 
POBOX  1146 
NORMAL  AL  35762 


ABERDEEN  PROVING  GROUND 

25  DIR  USARL 

AMSRL  WM  BD 
W  R  ANDERSON 
R  A  BEYER 
A  L  BRANT 
S  W  BUNTE 
C  F  CHABALOWSKI 
T  P  COFFEE 
J  COLBURN 
P  J  CONROY 
B  E  FORCH 
B  E  HOMAN 
S  L  HOWARD 
P  J  KASTE 
A  J  KOTLAR 
C  LEVERITT 
K  L  MCNESBY 
M  MCQUAID 
A  W  MIZIOLEK 
J  B  MORRIS 
J  A  NEWBERRY 
M  J  NUSCA 

R  A  PESCE-RODRIGUEZ 
G  P  REEVES 
B  M  RICE 
R  C  SAUSA 
A  W  WILLIAMS 
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